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MEASUREMENT  OF  GROUND  DIELECTRIC 
PROPERTIES  USING  WIDE-ANGLE 
REFLECTION  AND  REFRACTION 

Steven  A.  Arcone  and  Allan  ).  Delaney 


INTRODUCTION 

Knowledge  of  the  dielectric  properties  of  soils 
is  essential  to  the  interpretation  of  ground  radar 
profiling,  the  calculation  of  ground  reflection 
coefficients,  and  the  determination  of  soil  water 
content  by  airborne  remote  sensing  In  continu¬ 
ous  ground  radar  profiling  (Annan  and  Davis 
1976,  Davis  et  al  1976),  the  dielectric  properties 
determine  signal  velocity  and  are  therefore  es¬ 
sential  to  translating  echo  return  times  into  geo¬ 
logic  reflector  geometry  Radar  profiling  is  usu¬ 
ally  carried  out  at  a  fixed  separation  between 
the  transmitting  and  receiving  antennas,  while 
the  best  method  for  obtaining  dielectric  values 
is  to  vary  the  antenna  (geophone)  separation 

Ground  probing  radar  has  been  applied  most 
successfully  in  the  very  high  frequency  (VHF) 
band  of  30  to  300  MHz  At  high  frequency  (3  to 
30  MHz)  and  VHF,  horizontally  polarized  waves 
have  been  used  to  measure  ground  dielectric 
permittivity  and  conductivity  by  the  method 
known  as  radio  frequency  interferometry  (RFI) 
(Annan  1973)  In  this  method,  electric  properties 
are  calculated  from  the  interference  pattern  be¬ 
tween  continously  propagating  air  and  ground 
waves  The  method  has  been  used  at  frequencies 
as  low  as  2  MHz  (Rossiter  et  al.  1973)  and  as  high 
as  88  MHz  (Arcone  and  Delaney  1981). 

As  offshoot  of  the  RFI  method  that  employs 
short  pulses  was  discussed  by  Annan  and  Davis 
(1976)  In  this  technique,  known  as  WARR  (wide- 
angle  reflection  and  refraction),  the  time  delays 


between  air  ground  wave  pulses  are  measured  to 
determine  dielectric  properties.  Annan  and 
Davis  used  a  ground-probing  radar,  the  antennas 
of  which  were  continuously  separated  They  re¬ 
corded  their  results  on  strip  chart  forms  where 
the  horizontal  axis  was  antenna  separation  and 
the  vertical  axis  the  time  of  propagation  The 
slope  of  the  first  signal  received  after  the  air 
wave  then  gave  them  a  good  measure  of  the 
ground  surface  dielectric  constant. 

An  example  of  a  continuous  WARR  profile  is 
shown  in  Figure  1  The  dark  bands  represent  sig¬ 
nal  amplitude  and  the  thin,  light  bands  are  zero 
voltage  levels  in  the  return  signals  The  various 
sloping  profiles  correspond  to  a  variety  of  propa¬ 
gation  modes  (discussed  in  the  next  section).  In 
this  type  of  display,  information  that  could  be 
derived  from  the  waveforms,  such  as  dispersion, 
interference  and  phase  reversals,  is  not  avail¬ 
able 

In  this  report  the  ability  of  the  WARR  tech¬ 
nique  to  evaluate  dielectric  properties  is  thor¬ 
oughly  examined  by  displaying  the  individual 
waveforms  at  several  discrete  antenna  separa¬ 
tions  as  in  seismic  work  The  soundings  are  com¬ 
pared  with  continuous  radar  reflection  profiles 
at  fixed  antenna  separation  and  also  with  low 
frequency  resistivity  profiles  in  order  to  extract 
the  maximum  amount  of  dielectric  information 
For  the  WARR  soundings,  the  antennas  were  se¬ 
parated  in  1  0-m  intervals  to  14  or  1 5  m  to  inves¬ 
tigate  the  top  4  to  5  m  of  low  loss  sands  and  gra¬ 
vels  at  three  sites  in  interior  Alaska 


Figure  1.  Example  of  a  continuous  WARR  profile  I from 
Annan  and  Davis  1976).  The  horizontal  axis  is  antenna  sepa¬ 
ration  in  feet.  The  darkness  (paper  burn)  depends  on  signal 
intensity. 


THEORY  OF  GROUND  WAVE  PROPAGATION 
FROM  A  HORIZONTAL  ELECTRIC  DIPOLE 

The  theory  of  wave  propagation  in  the  steady 
state  from  a  horizontal  electric  dipole  lying  on 
or  near  the  earth's  surface  has  been  well  treated 
by  Kong  (1972),  Annan  (1973),  Annan  et  al  (1975) 
and  Chlamtac  and  Abramovici  (1981)  as  well  as 
by  many  others.  These  discussions  provide  a 
good  basis  for  understanding  pulse  propagation 
in  homogeneous  ground  because  pulses  may  be 
synthesized  from  an  infinite  number  of  plane 
wave  trains  In  layered  ground,  however,  this  ap¬ 
proach  will  not  apply  because,  with  continuous 
radiation,  waveguide  modes  exist  that  are  not 
analogous  to  the  events  observed  with  pulses. 

Tigure  2  (from  Annan  1973)  presents  the  four 
modes  that  are  excited  when  a  horizontal  elec¬ 
tric  dipole  is  placed  on  the  surface  of  a  homo¬ 
geneous  model  earth  The  modes  are  referenced 
to  a  right-hand  Cartesian  coordinate  system,  x,  y, 


z  In  the  figure,  wavefront  A  represents  a  spheri¬ 
cal  air  wave  above  the  surface,  and  wave  B  re¬ 
presents  a  spherical  wave  excited  in  the  ground. 
Waves  C  and  l;  propagate  above  and  below  the 
ground,  respectively,  in  order  to  maintain  the 
continuity  of  the  tangential  electric  and  magne¬ 
tic  fields  of  waves  B  and  A  across  the  air/ground 
interface.  Wave  C,  the  ground  surface  wave,  is 
termed  "inhomogeneous"  in  that  it  attenuates 
exponentially  with  height  above  the  surface  It 
propagates  with  the  same  horizontal  phase  velo¬ 
city  as  B.  Wave  D  is  called  a  head  wave  and  it 
propagates  with  the  same  horizontal  phase  velo¬ 
city  as  A  The  RFI  technique  is  based  on  the  in¬ 
terference  between  waves  C  and  A,  while  the 
WARR  technique  is  based  on  their  separation  in 
time 

The  velocity  of  wave  propagation  in  the 
ground  is  much  slower  than  in  air.  Consequently, 
in  order  for  A  to  match  D  at  the  air/ground  inter¬ 
face,  wave  D  must  propagate  at  an  angle  p  such 
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sumed.  The  loss  tangent  of  a  medium  is  defined 
as 

tan  d  =  (x"  +  'Uqwc0)/k'.  (4) 

In  air  the  phase  velocity  of  wave  A  is  simply  the 
phase  velocity  c  of  light  in  a  vacuum.  The  phase 
velocity  of  wave  C  is  w/fc,  which  for  tan  6  <  1, 
equals  cIs/HT,. 

The  addition  of  a  subsurface  interface  creates 
several  additional  events  as  illustrated  by  the 
rays  in  Figure  3.  The  layering  produces  a  wave¬ 
guide  effect  and  gives  a  1/W dependency  to  the 
amplitudes  of  the  guided,  subsurface  reflection, 
in  contrast  to  the  1/y2  dependency  for  the 
ground  surface  wave.  Of  course,  the  initial  am¬ 
plitudes  of  all  the  multiple  reflections  will  be  re¬ 
duced  from  that  of  the  ground  surface  wave.  If  a 
layer  is  sufficiently  thin  there  will  also  be  a  low 
frequency  filtering  effect  on  the  pulse,  but  this 
was  not  encountered  in  these  studies.  Also 
shown  in  Figure  3  is  the  possibility  that  a  re¬ 
fracted  wave  will  occur  if  x*2  <.  x*,. 


Figure  2.  Various  modes  radiated  in  the  broad¬ 
side  direction  from  a  horizontal  electric  dipole  ly¬ 
ing  on  the  surface  of  a  homogeneous  earth  ( from 
Annan  1973).  The  modes  are  described  in  the  text. 
The  angle  /I  =  sin  1  (k(/k,). 


Figure  3.  Possible  reflection  and  refraction 
events  in  a  two-layer  ground.  A  is  the  air  wave,  8, 
the  ground  surface  wave,  B,  and  B ,  reflections 
from  the  subsurface  interface  and  E  a  refracted 
wave  which  will  occur  if  ,x2:  ix,;. 


that  sin  (1  =  k0lk ,  where  k0  and  k,  are  the  propa¬ 
gation  functions  in  air  and  ground,  respectively. 
These  functions  are  given  by 

k0  =  (1) 

and 

k ,  =  0  +  iwMo/c  (2) 

where  x*  is  the  relative  complex  dielectric  per¬ 
mittivity,  which  may  be  expressed  as 

x*  =  x'  +  »x".  (3) 

x'  is  usually  termed  the  dielectric  constant  and 
x"  the  loss  factor  In  these  expressions  w  is  fre¬ 
quency  in  rad/s,  m»  =  4nx10  ’  H/m,  f„  = 
8  85x10  u  F/m,  t>  is  resistivity  in  ohm-m,  i  = 
and  the  time  dependence  e'"'’1  has  been  as¬ 


EQUIPMENT  AND  METHODS 

The  equipment  used  for  this  investigation  was 
a  Geophysical  Survey  Systems  Inc.  radar  system, 
featuring  a  model  4000  mainframe,  model  700P 
control  module  and  model  76  transmitting  and 
receiving  antennas.  The  antennas  (resistively 
loaded,  bowtie-type  dipoles)  were  in  two  separ¬ 
ate  housings.  The  transmitting  antenna  emits  a 
pulse  of  several  nanoseconds  duration,  and  the 
received  signals  are  sampled  and  converted  to 
an  audio  frequency  facsimile  for  display.  The 
time  scale  can  be  determined  from  a  supplied, 
calibrated  oscillator  or  by  physically  separating 
the  antennas  a  known  distance  and  calibrating 
the  time  delay,  as  was  done  for  this  investiga¬ 
tion. 

Figure  4a  shows  the  shape  of  the  pulse  emit¬ 
ted  when  the  antenna  radiates  into  space.  This 
pulse  shape  was  obtained  by  placing  the  anten¬ 
nas  over  a  perfectly  conducting  plane  and  separ¬ 
ating  them  several  free-space  wavelengths  of  the 
pulse  center  frequency  (1.0  m).  At  closer  dis¬ 
tances,  the  received  pulse  was  a  synthesis  of  di¬ 
rect,  indirect  (i.e.  reflected)  and  ground  waves 
By  about  3  m,  however,  the  direct  and  indirect 
waves  had  canceled  each  other  and  only  an  air 
surface  wave,  which  has  the  same  time  depend¬ 
ency  as  a  free  space  wave,  remained.  The  dura¬ 
tion  of  one  cycle  of  this  far-field  pulse  shape  is 
about  3.3  ns,  giving  a  center  frequency  to  the 
Fourier  spectrum  of  about  300  MHz. 

Figure  4b  shows  the  alteration  of  this  pulse 
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a.  Air  surface  wave  over  perfect  conductor. 


b.  Ground  surface  wave  in  moist  sandy  silt  for  which  x' 
s  70. 


Figure  4.  Far-field  radar  pulse  shapes. 


shape  when  the  antennas  are  coupled  to  ground. 
In  this  case  the  antennas  are  separated  over  a 
sandy  silt  with  a  dielectric  constant  of  about 
10.0.  The  frequency  content  of  this  ground  sur¬ 
face  wave  over  the  perfect  conductor  and  the 
relative  amplitudes  between  the  peaks  have 
changed. 

The  theoretical  difference  in  time  delay  be¬ 
tween  the  air  wave  signals  propagating  over  an¬ 
tenna  separations  of  1  and  14  m  was  used  to 
check  the  calibration  of  an  internally  generated 
time  scale  for  all  positions.  The  internal  time  cal¬ 
ibration  was  always  within  ±1.5  ns  of  the  theo¬ 
retical  difference  of  43.33  ns.  The  time  delays  of 
reflection  and  refraction  events  were  measured 
from  an  internally  recorded  "start  of  scan"  pulse 
which  initiated  every  scan  of  the  pulse  returns. 
The  returns  were  recorded  in  an  "A-scope"  (am¬ 
plitude  vs  time)  type  of  display. 

Two  different  antenna  positions,  both  produc¬ 
ing  horizontal  polarization  and  broadside  radia¬ 
tion.  were  used  for  the  WARR  profiling  at  each 
site.  One  was  the  normal  profiling  position  for 
which  the  antennas  were  designed.  This  position 
places  the  plane  of  the  "bowtie"  of  the  antenna 
parallel  with  the  ground  and  gives  the  best 
ground  coupling,  but  produces  an  air  wave  of 
very  small  amplitude  because  a  metal  shield  is 
attached  behind  the  antenna.  The  other  position 
simply  has  the  bowtie  turned  on  its  broadside 
edge.  This  position  gives  poor  ground  coupling 
but  the  air  wave  is  strongly  launched. 

To  verify  the  presence  of  any  suspected  layer¬ 
ing,  continuous  radar  profiles  were  performed 


along  each  sounding  line  of  the  WARR  measure¬ 
ments  using  a  different  antenna  set  (transmitter 
and  receiver  separated  by  30  cm)  housed  in  one 
unit.  Low  frequency  resistivity  measurements 
were  also  made  to  verify  the  dielectric  nature  of 
the  ground.  The  resistivity  measurements  were 
performed  with  a  magnetic  induction  device 
(Geonics  EM-31),  the  theory  and  operation  of 
which  may  be  found  in  Arcone  et  al.  (1979)  and 
Arcone  and  Delaney  (1981).  The  device  uses  two 
loop  antennas  separated  by  3.66  m  and  can  em¬ 
ploy  them  in  two  different  orientations  that  al¬ 
low  approximate  penetration  depths  of  3  and  6 
m.  Increases  in  d.c.  resistivity  were  interpreted 
as  indicative  of  decreasing  water  content  and 
vice  versa.  This  interpretation  implies  a  corre¬ 
sponding  decrease  in  dielectric  constant  at  VHF 
_with  increasing  resistivity. 


RESULTS 
Site  1 

This  site  is  part  of  a  glacial  outwash  plain  of 
sand,  silt  and  gravel  located  at  Ft.  Greely, 
Alaska,  and  now  used  infrequently  as  an  airfield. 
The  area  is  flat  and  the  organic  cover  has  been 
cleared.  Immediately  below  the  surface  are  nu¬ 
merous  fossil  ice  wedge  casts  filled  with  silt.  The 
horizontal  beds  of  stratified  sand  and  gravel 
about  the  casts  are  generally  distorted  towards 
the  vertical.  Church  et  al.  (1965)  have  docu¬ 
mented  many  cross  sections  of  these  casts.  A 
perched  water  table  exists  at  about  10  m. 
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a.  Low  frequency  magnetic  induction  resistivity  profile 
for  penetration  depths  of  3  and  6  m. 


Distance  (m) 


fa.  Radar  profile. 


Figure  5.  Continuous  radar  reflection  and  low  frequency 
magnetic  induction  resistivity  profiles  along  an  18-m 
traverse  at  Ft.  Creely. 


Figure  5  shows  the  continuous  radar  and  resis¬ 
tivity  profiles  along  an  18-m  segment  of  the  air¬ 
field.  The  resistivity  profile  shows  a  surprising 
lack  of  variation  both  vertically  and  horizontal¬ 
ly,  but  the  dielectric  contrast  between  the  silt  in 
the  casts  and  the  sand  and  gravel  causes  the 
complexity  of  the  radar  profile.  These  strong  re¬ 
turns  mainly  show  multiple  hyperbolas  caused 
by  scattering  from  the  localized  casts.  Further 
discussion  of  this  and  continuous  radar  profiles 
from  this  site  will  be  presented  in  Arcone  et  al. 
(1982). 

Figure  6  shows  the  individual  WARR  traces  ta¬ 
ken  across  the  profile.  The  transmitting  antenna 
was  situated  at  the  0-m  station  of  Figure  5  and 
the  receiving  antenna  at  stations  1-15.  Times  at 
which  the  events  of  strongest  amplitude  align 
from  trace  to  trace  are  connected  by  straight 
lines  and  replotted  in  Figure  7  The  slopes  of 
these  plots  allow  the  signal  velocity  to  be  calcu¬ 
lated  and.  hence,  also  the  corresponding  dielec¬ 
tric  constant.  The  plot  at  the  top  of  Figure  7  is 
for  the  air  wave,  and  the  data  points  for  this  plot 


were  taken  from  a  separate  WARR  spread  with 
the  antennas  oriented  on  edge  for  better  air 
wave  coupling.  The  intersection  of  the  air  wave 
plot  with  the  vertical  axis  gives  the  correct  zero 
time  reference 

The  slope  of  each  plot  gives  a  dielectric  con¬ 
stant  of  about  5.5.  Since  x'  for  dry  sand  or  gravel 
can  be  as  low  as  2.6,  some  moisture  was  in  the 
ground,  as  also  shown  by  the  low  resistivity  val¬ 
ues.  The  later  events  shown  in  Figure  7  are  be¬ 
lieved  to  have  scattered  from  a  silt-filled  cast  in 
the  vicinity  of  the  transmitter  (0  m)  and  are  re¬ 
vealed  by  the  wide  hyperbola  starting  downward 
at  about  0  and  2  m  at  20  ns  in  Figure  5.  If  these 
plots  were  of  multiple  reflections  from  a  layer 
interface,  then  their  slopes  would  have  appeared 
to  be  merging  asymptotically. 

These  WARR  dielectric  results  give  the  same 
implication  of  homogeneity  as  the  resistivity 
profile,  probably  because  both  were  integrating 
electrical  responses  that  occurred  over  horizon¬ 
tal  and  vertical  distances  that  were  greater  than 
those  of  the  casts  (approximately  1-2  m).  The 
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Figure  6  WARR  sounding  at  Ft. 
Greely.  The  first  pulse  at  0  ns  is  arti¬ 
ficially  generated  to  start  the  scan 
and  serves  as  a  time  reference. 


Figure  7  Time-distance  plots  taken  from 
some  of  the  events  observed  in  f  igure  b 


continuous  radar  profile,  however,  does  not  inte¬ 
grate  over  large  distances  so  that  greater  detail 
was  revealed 

Site  2 

This  site  is  located  at  Ft  Wainwright  in  Fair¬ 
banks.  Alaska,  and  is  part  of  the  Tanana  River 
floodplain.  The  site  material  is  mainlv  silt,  sand 
and  gravel,  with  most  of  the  silt  having  been 
bladed  off  to  form  a  road.  The  silt  is  5.1  m  above 
the  Chena  River  which  was  within  TOO  m.  Thr  e- 
fore  the  water  table  b  at  a  depth  of  approximate¬ 
ly  5.1  m 

Figure  8  shows  the  continuous  radar  and  re¬ 
sistivity  profiles  along  a  14-m  section  of  the 
road  The  second  layer  occurs  at  a  delay  of  35  to 
40  ns  and  is  probably  the  base  of  the  frozen  ma¬ 
terial  The  third  layer  seems  to  occur  at  a  delay 
of  about  65  ns  and  presents  a  more  diffused  re¬ 
turn.  This  event  is  assumed  to  be  the  water  table. 
The  resistivity  profile  gives  much  higher  values 
than  at  site  1  and  hows  resistivity  to  be  increas¬ 
ing  with  depth  This  increase  may  be  due  to 
greater  amounts  of  coarse-grained  material  at 
depth  or  to  a  decrease  in  moisture  content  down 
to  the  water  table,  as  might  be  caused  by  up¬ 
ward  moisture  migration  during  freezing.  Unfro¬ 
zen 'adsorbed  water  is  presumed  to  cause  the 
lowest  resistivities  near  the  surface. 
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figure  8  Continuous  radar  reflection  and  low  fre¬ 
quency  magnetic  induction  resistivity  profiles  along 
a  N-m  traverse  at  ft.  Wainwnght  ibladed  road  site/. 


Figures  9  and  10  show  the  individual  WARR 
sounding  traces  (transmitter  at  0  m)  and  the  cor¬ 
responding  time-distance  plots  for  several  of  the 
events.  The  air  wave  shown  here  and  the  ground 
surface  wave  do  not  intercept  the  time  axis  at 
the  exact  same  point  because  the  antennas  were 
in  the  normal  profiling  position,  which  does  not 
allow  good  air-wave  coupling.  The  other  plots 
show  a  very  characteristic  two-layer  behavior 
with  the  slopes  all  merging  at  large  distances  As 
it  is  impossible  to  pick  out  the  leading  edge  of 
each  multiple  event,  the  depth  of  the  first  layer 
interface  is  estimated  to  be  about  2.5  m  and  the 
dielectric  constant  of  this  first  layer  is  calcu¬ 
lated  from  the  slopes  of  the  plots  to  be  about 
4.8 

Since  the  water  table  is  at  5.1-m  depth,  the  se¬ 
cond  layer  is  about  2.6-m  thick.  The  second  layer 
return  is  about  30  ns  behind  the  first  layer  return 
in  Figure  8  so  that  the  dielectric  constant  of  the 
second  layer  is  calculated  at  about  3.0  using  the 
formula  x'  =  ( ctl2d)‘  where  d  is  layer  depth,  c  is 
free  space  velocity  and  t  is  round-trip  time  of 
propagation  within  the  layer.  This  value  is  con¬ 


sistent  with  the  above  hypothesis  of  a  drier  ma¬ 
terial  as  deduced  from  the  resistivity  data  This 
value  also  implies  that  a  refracted  wave  may  oc¬ 
cur  along  the  interface  between  layers  1  and  2, 
as  the  second  layer  is  calculated  to  be  of  a  high¬ 
er  velocity  medium  than  the  first 

This  analysis  provides  the  subsurface  model 
for  site  2  shown  in  Figure  11  Also  shown  in  Fi¬ 
gure  11  is  a  ray  path  for  a  refracted  ray  which 
should  exist  for  this  model.  The  angle  of  inci¬ 
dence  upon  the  first  subsurface  interface  for  the 
refracted  ray  is  sin'1'  \/ k',Ik' ,  =  52.2°.  The  first 
antenna  station  to  receive  this  refracted  ray 
would  be  at  7.0  m,  as  the  refraction  first  reaches 
the  surface  at  the  theoretical  distance  of  about 
6.5  m  At  the  7-m  station,  the  refracted  wave 
should  be  about  15  ns  behind  and  in  phase  with 
the  ground  surface  wave  An  event  thought  to  be 
this  refraction  occurs  at  the  7-m  station  in  Figure 
9  The  exact  phase  agreement  between  this 
event  and  the  ground  surface  wave  is  not  evi¬ 
dent  since  the  pulses  overlap.  The  ground  sur¬ 
face  and  refracted  waves  should  then  cross  over 
(in  phase  with  no  time  displacement)  at  the  15-m 
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Figure  9.  Radar  WARR  sounding,  Ft.  Wainwright  (bladed  road 
site). 
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Figure  W.  Time-distance  plots  taken  from 
events  observed  in  Figure  8. 


station.  At  the  14-m  station  the  theoretical  time 
delay  difference  is  4  ns  (or  just  slightly  less  than 
one-half  the  period  of  the  pulse)  and,  in  fact,  ser¬ 
ious  distortion  of  the  leading  edge  of  the  ground 
surface  wave  occurs  here. 

Site  3 

This  site,  also  located  at  Ft  Wainwright,  is  in 
the  same  area  as  site  2  and  of  the  same  kind  of 
gravelly  sand.  The  area  has  been  used  as  a  gravel 
pit  No  visible  evidence  of  the  water  table  was 


present  in  nearby  pits  but  it  is  known  to  be  more 
than  5  m  deep. 

Figure  12  shows  the  continuous  radar  and  re¬ 
sistivity  profiles  along  a  14-m  section  of  the  area. 
The  radar  profile  shows  a  possible  two-layer 
structure  between  0  and  8  m  and  then  some  un¬ 
interpretable  returns  between  8  and  14  m  The 
resistivity  profiles  show  that  the  near-surface  re¬ 
sistivity  more  than  doubles  across  the  section 
while  the  deeper  values  generally  stay  between 
600  and  900  ohm-m. 

The  individual  WARR  sounding  traces  in  Fi¬ 
gure  1 3  also  show  a  severe  disturbance  in  electri¬ 
cal  structure  beginning  at  7  m.  More  data  could 
have  been  collected  to  resolve  this  structure, 
but  the  primary  purpose  was  to  measure  dielec¬ 
tric  properties  In  Figure  14  the  earliest  returns 
give  a  dielectric  constant  of  2  9,  which  implies  a 
very  dry  surface  layer  of  sand  or  gravel.  The  low 
value  is  most  usually  associated  with  resistivity 
values  at  10*  ohm-m  or  above  A  secondary 
wavelet  gives  a  dielectric  value  of  7  4  out  to  a 
distance  of  9  m  before  this  return  is  interfered 
with  by  the  returns  from  the  inhomogeneity  The 
wavelet  traveling  at  k'  -  2.9  is  the  ground  sur¬ 
face  wave,  whereas  the  wavelet  traveling  at  *' 

7  4  is  probably  the  part  of  the  ground  surface 
wave  wavefront  that  has  been  continuously  re¬ 
fracted  downwards  and  reflected  from  the  inter¬ 
face  shown  at  about  35  ns  in  Figure  12  Within 
this  layer  the  dielectric  constant  probably  in¬ 
creases  continuously  with  depth  If  7.4  is  about 
the  average  dielectric  constant  for  the  layer, 
then  a  time  delay  of  about  35  ns  gives  a  layer 
depth  of  about  19  m 
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Figure  11.  Subsurface  model  for  site  2  based  on  the  WARR  and 
continuous  radar  profiling  data.  The  surface  elevation  of  the 
site  above  the  nearby  Chena  River  was  used  for  the  water  table 
depth  of  5.1  m,  the  WARR  data  were  used  for  calculating  k ,  = 
4.8  and  the  first  layer  depth  of  2.5  m,  and  the  continuous  radar 
reflection  profile  was  used  to  calculate  =  Ft) 
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Figure  12.  Continuous  radar  reflection  and  low  fre¬ 
quency  magnetic  induction  resistivity  profiles  along 
14-m  traverse  at  site  3  (Ft.  Wainwright  gravel  pit  site l 


SUMMARY  AND 
CONCLUDING  REMARKS 

These  investigations  have  shown  that  WARR 
soundings  are  an  effective  aid  to  radar  profiling 
because  of  the  dielectric  information  they  pro¬ 
vide  At  site  1  the  *'  value  of  5  5  obtained  for  the 


sand  and  gravels  was  shown  to  be  an  effective 
bulk  value  over  the  few  meters  of  depth  of  the 
silt  ice  wedge  casts.  At  site  2  the  combination  of 
WARR,  continuous  profiling  and  knowledge  of 
the  water  table  allowed  the  dielectric  permittivi¬ 
ty  of  two  separate  layers  to  be  evaluated.  At  site 
3,  the  WARR  soundings  revealed  that  the  bulk 
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Figure  14.  Time-distance  plots  taken  from 
events  observed  in  Figure  12. 


permittivity  value  of  the  upper  few  meters  of 
ground  was  not  the  same  as  the  surface  value  so 
that  permittivity  was  increasing  with  depth. 

At  all  three  sites  the  ground  surface  was  flat 
so  that  the  antennas  were  in  close  ground  con¬ 
tact  and  no  elevation  corrections  had  to  be  ap¬ 
plied  to  the  traces.  If  the  antennas  had  been  ele¬ 
vated  by  any  surface  roughness  or  by  a  snow- 
pack,  the  ground  surface  wave  events  would 
have  been  seriously  attenuated.  This  is  because 
the  amplitude  dependency  on  height  above 
ground  for  a  horizontally  polarized  ground  sur¬ 
face  wave  is  proportional  to  exp  [-(kj  -  fcj)'1  z) 
where  z  is  height.  At  the  value  x'  =  5,  this  attenu¬ 
ation  rate  works  out  to  be  about  0  25  dB/cm  at  a 
frequency  of  80  MHz. 

The  dielectric  properties  of  soils  show  little 
dispersion  for  small  water  contents  (i.e.  less  than 
5%  by  volume),  as  was  shown  by  the  uniformity 
of  shapes  seen  from  trace  to  trace  in  Figures  6,  9 
and  13  for  the  ground  surface  waves  Although 
the  relative  amplitudes  of  the  peaks  changed 
slightly,  the  period  of  oscillation  remained  fairly 
constant  throughout  the  sounding  Dispersion  is 
associated  with  high  loss  tangents  The  highest 
possible  loss  tangent  that  may  have  been  en¬ 
countered  would  have  been  at  the  Ft  Creely  air¬ 
strip  where,  for  a  d  c.  resistivity  of  about  300 
ohm-m  and  a  wavelet  center  frequency  of  80 
MHz,  the  loss  tangent  (eq  4)  would  have  been 
moderate  at  about  0.13  However,  in  this  case 
the  ground  surface  waveforms  were  consistently 
uniform  with  little  dispersion. 
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The  main  advantage  of  continuous  WARR  pro¬ 
filing  is  that  reference  points  within  the  subsur¬ 
face  returns  may  be  easily  traced  to  allow  the 
calculation  of  dielectric  constants  and  layer 
depths  When  used  in  the  discrete  mode,  the 
waveform  display  allows  dispersion,  interfer¬ 
ence  and  phase  changes  to  be  identified.  In  addi¬ 
tion,  good  antenna-ground  coupling  need  only 
be  maintained  at  discrete  stations  Finally  it 
should  be  mentioned  that,  in  the  discrete  mode, 
signal  enhancement  techniques  such  as  stacking 
and  deconvolution  can  be  applied  most  easily. 
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